Singularity free stars in (2+1) dimensions 
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We present some new types of non-singular model for anisotropic stars with constant A and 
variable A based on the Krori and Barua (KB) metric in (2+ 1) dimensions. The solutions obtained 
here satisfy all the regularity conditions and its simple analytical form helps us to study the various 
physical properties of the configuration. 



PACS numbers: 04.40.Nr, 04.20.Jb, 04.20.Dw 
I. INTRODUCTION 

The study of (2+1) dimensional gravity has become a 
subject of considerable interest. Newtonian theory can 
not be obtained as a limit of Einstein's theory in (2+1) 
dimensional spacetime. In this case gravity is localized 
that means there is no propagation of gravity outside 
the sources. Also, it is argued that (2+1) dimensional 
gravity provides some new features towards a better un- 
derstanding of the physically relevant (3+1) dimensional 
gravity [1-4]. Most of the studies on this respect are black 
hole spacetimes, star or cosmological models [5]. How- 
ever, in recent years various studies have been done. Ra- 
haman et al[6] have proposed a new model of a gravastar 
in (2+1) anti-de Sitter space-time. They have also gen- 
eralized their earlier work on gravastar in (2 + 1) dimen- 
sional anti-de Sitter space-time to (2 + 1) dimensional 
solution of charged gravastar [7]. Some authors [8] have 
also discussed wormhole solutions in (2+1) dimensional 
spacetime. 

In recent past, Krori and Barua (KB) [9] constructed 
static, spherically symmetric solutions based on a par- 
ticular choice of the metric components 500 and gn in 
curvature coordinates. Recently, KB's approach was 
adopted by various authors for constructing star mod- 
els [10-15]. These studies are confined within (3+1) di- 
mensional spacetimes. Therefore, it will be interesting 
to search whether nonsingular solutions of KB type will 
be existed in (2+1) dimensional fluid sphere. We are 
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looking forward to get some extra features as our (2+1) 
dimensional models include an additional parameter, cos- 
mological constant. In 1999, Lubo et al [2] has discussed 
regular (2 + 1) spherically symmetric solutions, however, 
there approach were different. Anisotropy was used in 
the compact star configuration to allow some interesting 
studies and we mention some papers that yield meaning- 
ful results [16]- [19]. 



In this investigation, we explore the possibility of ap- 
plying the Krori-Barua [9] metric to describe the interior 
spacetime of a star in (2+1) dimension. Nowadays, it is 
known that the dark energy represents 73% of the whole 
mass-energy of our Universe. This conclusion is given 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
that also indicates that the dark-energy is causing a 
speeding of the expansion of the rate of the universe. 
The subject of compact stars is of actuality and under 
study in the last decades. One of the possibility for the 
formation of compact anisotropic stars is to use the cos- 
mological constant. For this reason, we have considered 
cosmological constant in our model. We have discussed 
two models, one with constant A and the other with vari- 
able A. 



The structure of our work is as follows: In section II, 
the non-singular model for anisotropic stars with con- 
stant A based on the Krori and Barua (KB) metric in 
(2 + 1) dimensions is developed. In section III, we have 
discussed some physical features of the model. In section 
IV, we have presented the model with variable A. In Sec- 
tion V, we have analyzed some physical properties of the 
model given in section IV. Finally, in section VI, we have 
made a conclusion about our work. 
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II. NON-SINGULAR MODEL FOR 
ANISOTROPIC STARS WITH CONSTANT A 

Let us assume that the interior space-time of a star is 
described by the KB metric 



(1) 



with 2/i(r) = Ar"^ and 2v{r) = Br"^ + C where A, B and 
C are arbitrary constants which will be determined on 
the ground of various physical requirements. The energy 
momentum tensor in the interior of the anisotropic star 
is assumed in the following standard form 

Tij = diag{p, -pr, -pt), 

where p, Pr and pt correspond to the energy density, nor- 
mal pressure and transverse pressure respectively. 

Therefore, the Einstein field equations for the metric 
(1) with constant A can be written as (assuming natural 
units G = c = 1) 



27r/9-|- A 



27rpr — A 



p, e 



-2ij, 



V e 



-2n 



2-Kpt-K = e-^'' {y" + y'^ - y' ijl') . 



(2) 

(3) 
(4) 



Now, plugging the metric (1) in equations (2) - (4), we 
get the following expressions of energy density p, normal 
pressure Pr , tangential pressure pt as 



A. Matching Conditions 

The exterior (p = p = 0) solution corresponds to a 
static, BTZ black hole is written in the following form as 

ds2 = _ (_Mo - Ar^) de + (-Mo - Ar^)"' dr^+r^d0\ 

(10) 

Here we match our interior metric to the exterior BTZ 
metric and as a consequence, we get 



A = -^\niK^R^ - Mo), 



C-- 



B 



K'^B? - Mo ' 



K'^E? - Mo 



ln{K^R^ - Mo). 



(11) 
(12) 

(13) 



B. Regularity at the centre 

In this analysis, we consider anti-de Sitter space-time 
and to ensure that cosmological constant is always nega- 
tive, we write A = —K"^. Since the radial EOS is always 
less than unity, therefore equation (8) at once indicates 
that the cosmological constant should be negative. 

Now, we find the central density and central pressures 
( radial and transverse ) as 



Po = p{r = 0) = ^{A + K^) 



(14) 



1 

" = 2^ 



Ae-^^ - A 
Be-^^^ + A 



Pt 



I 

2^ 



e-^'- (5^2 +B- ABr'') + A 



(5) 
(6) 
(7) 



p,{r = 0)=pt{r = 0)=po = — {B-K^) (15) 
Prom the equations (5) and (6) , we find 



dp 



A^ 



, = re 

dr IT 



<o, 



Using equations (5) - (7), the equations of state (EOS) 
parameters corresponding to radial and transverse direc- 
tions are written as 



Wr(r) = 



Be 



-At'- 



■A 



A^-Av^ _ A 



(8) 



and 



dpr 
dr 



AB 



re-^'^ <0. 



which gives a meaningful result that density and pressure 
are decreasing function of r. The above equations imply 
at r = 0, 



wt(r) 



-Ar^ (B^r^ + B - ABr^) + A 
Ae-^"-^ - A ■ 



(9) 



III. PHYSICAL ANALYSIS 

In this section we will discuss the following features of 
our model: 



and 



dp 
dr 



0, 



fp 
dr^ 



d'^Pr 

dr'^ 



dpr 

dr 



0, 



<0, 



AB 



< 0. 



(16) 
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which support maximahty of central density and radial 
central pressure. 

For our model, the measure of anisotropy, A — Pt—Pr, 
is given by 



[A/— 2.25soIar mass, R — \Okm ] 



A 



2tt 



(17) 



The anisotropy, as expected, vanishes at the centre. 

The energy density and the two pressures are contin- 
uous function of radial coordinate r that means they are 
well behaved in the interior of the stellar configuration. 
The radius of star is obtained by letting pr{r = R) = , 
which gives 



R 




(18) 



0.0085 



0.0080 



0.0075 



0.0070 




A = -0.035 A = -0.038 A = -0.04 



We have considered the data from a 3 spatial dimen- 
sional object in order to fix the constants of the model. 
Although it is not very clear the physical meaning by 
considering stelar objects of 2 spatial dimensions, how- 
ever, if an observer is sitting in the plane, 9 = constant 
, then he sees all characteristics as a (2+1) dimensional 
picture. So apparently, one can use all the data which 
are more or less the same for both spacetimes. 

Therefore we have used the data from X ray buster 
4U 1820-30 to calculate the corresponding constants. It 
is known that the mass of X ray buster 4U 1820-30 is 
2.25 Mq and radius 10 km. To understand the physi- 
cal behavior of the solutions of our model, we have as- 
sumed certain values of A. As we considered the anti 
de-Sitter spacetime, we have used negative values, say, 
A = -.035, -.038, -0.04 ( see reference [20] ). 

Using the data from X ray buster 4U 1820-30, we have 
obtained the values of the constants A and B via equa- 
tions (11) and (12) for different values of A in units of 
km^^ ( see table 1). 

Plugging in G and c in the relevant equations, we have 
calculated the central density po = 9.49 x lO^^gm cm~^, 
surface density p/j = 9.23 x 10^^ gm cm~^, central pres- 
sure Pr{r = 0) = pt(r = 0) = 3.28 x 10"^^ dyne cm^^ for 
A = -.04. 

We draw the figures for the density variation and pres- 
sures at the stellar interior of X ray buster 4U 1820-30 ( 
see figures 1 and 2 ). The figure 3 indicates that A > 

1. e. the 'anisotropy' is directed outward. Thus our model 
exerts a repulsive 'anisotropic' force (A > 0) which allows 
the construction of more massive distributions. In figure 

2, one can note that transverse pressure is increasing to- 
wards the surface. These types of stars are composed 
of unknown matters whose density is above the normal 



nuclear density, p„ ~ 4.6 x 10"'^* 



gm cm 



As a result. 



FIG. 1: Density ( along vertical axis ) variation at the stellar 
interior of X ray buster 4U 1820-30 of mass 2.25 Mq and 
radius 10 km for different values of A. 



0.12- 



0.10 



0.08 



[M ^2.25 solar mass, \ Okm] 
/' 



Pi 0.06 



0.04 



0.02 




these highly compact stars are anisotropic in nature and 
peculiar phenomena may be happened there. 



FIG. 2: Radial and transverse pressures variation at the stel- 
lar interior X ray buster 4U 1820-30 of mass 2.25 Mq and 
radius 10 km for different values of A. Solid lines and chain 
lines indicate radial and transverse pressures respectively. 



C. TOV Equation 

The generalized TOV equation for an anisotropic fiuid 
distribution is given by 

^ '^^^"'^^^ ^''^ + l^P^' P*) = 0- (19) 
According to Ponce de Leon [21] suggestion, one can 
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TABLE I: Values of the constants A, B for X ray buster 4U 1820-30 for different values of A. 



A 


M 


R (km) 


A 


B 


-0.035 


2.25 Mq 


10 


0.017078 


0.19310 


-0.038 


2.25 Mq 


10 


0.007313 


0.07896 


-0.04 


2.25 Mq 


10 


0.003834 


0.05869 



[M= 2.25 solar mass, R = 10 km 




FIG. 3: The anisotropic behaviour at the stellar interior X 
ray buster 4U 1820-30 of mass 2.25 AIq and radius 10 km for 
different values of A. 




FIG. 4: Three different forces acting on fluid elements of X 
ray buster 4U 1820-30 of mass 2.25 A/q and radius 10 km in 



static equilibrium is shown against r with A 



.04. 



Fg - -Br{p + pr), (23) 



rewrite the above TOV equation as d f A 



dr \ 2tt 

Majp + Pr) li^ d f A \ 1 _ ^^g-Ar^ (24) 

Fa = -{pt~Pr). (25) 

r 

where Mq = Mcij) is the gravitational mass inside a 
sphere of radius r and is given 

by The profiles of Fg, Fh and Fa for our chosen source 

are shown in Fig. 4. The figure provides the information 
about the static equilibrium due to the combined effect of 
Mc{r) = r^e^~v' , (21) pressure anisotropy, gravitational and hydrostatic forces. 



[ This can be derived from the Tolman-Whittaker formula 
using Einstein field equations.] 

This modified form of TOV equation delineates the 
equilibrium condition for the X ray buster 4U 1820-30 
subject to the gravitational and hydrostatic plus another 
force due to the anisotropic nature of the stellar object. 
Now, we write the above equation as 

Fg + Fh + Fa = 0, (22) 



D. Maximum mass-radius relation 

We calculate the mass m(r) within a radial distance r 

as 

m(r) = / 27rpfdf = -[1 + K^r^ - e'^''^]. (26) 
Jo 2 

In Fig. 5, we have plotted this mass to radius relation. 
One can note that the upper bound on the mass in our 
model can be written as 



where 
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[M = 2.25 solar mass, R= 10 km] 



m{r) 




FIG. 5: Variation of mass function at stellar interior for dif- 
ferent values of A. 



which implies 



m(r) 



The plot against r (see Fig. 6) indicates that the ra- 
tio '"^''^ is an increasing function of the radial parameter. 
It is interesting to note that the constraint on the max- 
imum allowed mass-radius ratio in our case falls within 
the limit to the (3-f 1) dimensional case of isotropic fluid 

sphere i.e., ( ^^^^) — 0.216 < | ( here we have used 

V / max 

A = -.04 ). 



M 1 + K'^R^ - e 
'r - 2R 



-AR-' 



(28) 



E. Compactness and redshift 

From the above mass function (26), we obtain the com- 
pactness of the star as 



u = — = — 1 + K^r^ -e 
r 2r \ 



(29) 



and correspondingly the surface redshift (Zg) is given by 
= (1-2m)-^-1, (30) 

where 



m{r) 



0.20 

0.18 

0.16- 

0.14 

0.12 



O.IOH 

0.08 

0.06 

0.04H 

0.02 



[M = 2.25 solar mass, R = 10km] 



V 



// 







4 6 
r{km 1 



10 



■ A = -0.035 A = -0.038 A = -0.04 



FIG. 6: Variation of ^^^^^ is shown against r for different 
values of A. 



[M = 2.25solar mass, R - lOkm 



0.20- 



0.15- 




FIG. 7: The variation of redshift function Zs is shown against 
r for different values of A. 



1 - i (l + X V - 



- 1. 



(31) 



Thus, the maximum surface redshift of our (2+1) dimen- 
sional star of radius 10 km can be found as Zg — 0.328 ( 
here we have used A = — .04 ). 

The figure 7 indicates the variation of redshift function 



Zs against r for different values of A. 
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IV. VARIABLE A 



and 



Now, we consider the model taking the cosmological 
constant as radial dependence i.e. A = A(r) = A^. 
(say). 

To get the physically acceptable stellar models, we as- 
sume that the radial pressure of the compact star is pro- 
portional to the matter density i.e. 



Pr = mp, m > 0, 



(32) 



where m is the equation of state parameter. 

Now, using the above equation of state and equations 
(1) - (4), we get the energy density p, normal pressure 
Pr, tangential pressure pt and cosmological parameter A^, 
respectively as 



P 

Pr 

Pt 



{A + B) _^^2 



27r(m + 1) 
m{A + B) _^^2 
27r(TO + 1)^ 



>0, 



2tt 

,-Ar' 



[Br^{B -A) + 



m(A + B) 
1 + m ^ 



m + 



- [mA - B\ . 



(33) 

(34) 

(35) 
(36) 



Also, the equation of state (EOS) parameters corre- 
sponding to normal and transverse directions can be writ- 
ten as 



Ulr{T) = m, 



ujt{r) =m + 



{m + l)(B -A)Br'^ 
{A + B) 



(37) 



(38) 



V. PHYSICAL ANALYSIS 

In this section we will discuss the following features of 
our model: 

One can see from equations (33) and (34) that 



and 



dp 
dr 



dpr 
dr 



A{A + B) ^^_^^. 
7r(TO + 1) 



<0, 



m{A + B)A^^_^^2 



7r(m + 1) 
Also, at r = 0, our model provides 



< 0. 



dp 
dr 



= 0, 



dpr 

dr 



= 0, 



(fp _ A{A + B) 
dr'^ n{m + 1) 



<0, 



d^Pr 

dr"^ 



m,A(A + B) 
Tr{m + 1) 



<0. 



which indicate maximality of central density and central 
pressure. 

The central density and radial central pressures are given 
by 



Po = p{r = 0) 



(A + B) 
2TT{m+l) 



Pr{r = 0) = pt{r = 0) = Po = 



m{A + B) 
27r(m -I- 1) 



(39) 



(40) 



It is obvious that the radial pressure vanishes at the 
surface i.e. at r=R, pr{r = R) = where R is the radius 
of the star. Equation (34) implies either m = or A 
= - B. But, A = - B implies p = 0, which is not possi- 
ble, therefore, we should take m=0. In other words, in 
(2+1) dimensional KB spacetime with variable A, only 
dust anisotropic model exists. 

The measure of anisotropy, pt—Pr, in this model 
is obtained as 



^_ B{B-A) ^,^_^^. 
27r 



(41) 



At the centre A = i.e. anisotropy dies out as expected. 

We match the interior metric to the exterior BTZ ma- 
tric and get fortunately the same values of constants A, 
B and C given in equations (11)-(13). 

As before, we have considered X ray buster 4U 1820- 
30 and have chosen the values of the parameters A 
and B as A=0.003834 and B = 0.05869 to obtain the 
central density and surface density which are given by 
Po = 1-36 X lO^^gm cm,-^ and pr = 9.234 x lO^^gm cm"^ 
respectively. 

In this case the TOY equation is written as 



Fg + Fh + Fa= 0, 



(42) 



where 



Fa 

Fa 



-Br{p), 

d f kr 

~dr V 2^ 



(43) 



We note that here the variable A contributes to the hy- 
drostatic force and transverse pressure provide the effect 
of pressure anisotropy. 

The profiles of Fg, F^ and Fa for our chosen source 
are shown in Fig. 8. The figure provides the information 
of the static equilibrium due to the combined effect of 
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[M = 2.25solar mass, R=Vikm,A = 0.003834 B = 0.05869] 




FIG. 8: Three different forces acting on the fluid elements of 
X ray buster 4U 1820-30 of mass 2.25 Mq and radius 10 km 
in static equiUbrium is shown against r. 



pressure anisotropy, gravitational and hydrostatic forces. 

The effective mass m{r) within a radial distance r is 
defined as 



me//(r) = 
Hence, we get 

me//(r) 



27r 



ZTT 



rdf 



(44) 



In Fig. 9, we plot this mass to radius relation. One can 
note that the upper bound on the mass in our model can 
be written as 



2me//(r; 
which implies 



1 



Meff ^ I 



^AR^ 



R 



2R 



(45) 



(46) 



The plot "^^ff'^^^ against r (see Fig. 10) indicates that 

the ratio ^''^ is an increasing function of the ra- 
dial parameter. Again, we see that maximum allowed 
mass-radius ratio in our case falls within the limit to 
the (3-|-l) dimensional case of isotropic fluid sphere i.e., 

( "^y^'"^ ) = 0.01598 < I ( here we have used 

V / max 

A=0.003834 and B = 0.05869 ). 

The compactness of the star is given as 



m(r) 



1 



1 - e 



(47) 



r 2r 

and correspondingly the surface redshift (Zs) is given by 
Zs^{l-2u)-i -1, (48) 



[M= l.lisolar mass, R = 10 km, A = 0.003834] 
0.15- 



0.10 



0.05 




FIG. 9: Variation of mass function against r at stellar interior. 



M = 2.25solar mass, R = 10 km, ^ =0.003834] 




FIG. 10: Variation of — is shown against r at stellar inte- 
rior. 



where 



1--(1 

r 



-Ar' 



(49) 



Thus, the maximum surface redshift of our (2+1) dimen- 
sional star of radius 10 km can be found as Zs — 0.0164 
( here we have used A=0.003834 and B = 0.05869 ). The 
figure 11 indicates the variation of redshift function Zs 
against r. 
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[M = 2.2isolar mass, R = 10 km, ^ = 0.003834] 




2 4 6 8 10 



for stars with constant A and variable A. In this phe- 
nomenological model, we have employed physical data of 
the X ray buster 4U 1820-30 in our models. It is found 
that the central density po — 9.49 x lO^^^m cm^'^ and 
surface density pa — 9.23 x IQ^^gm cm~^ for the first 
case which are beyond the normal nuclear density. For 
the later case, densities are more than the former one. 
The models are attainable in static equilibrium condi- 
tions due to the combined effect of pressure anisotropy, 
gravitational and hydrostatic forces. One can note that 
the maximum allowed mass-radius ratio in our cases falls 
within the limit to the (3-1-1) dimensional case of isotropic 

fluid sphere i.e., ( '""y j < |. We have seen that 

V / max 

in (2-1-1) dimensional KB spacetime with variable A, only 
dust anisotropic model exists. This feature is absent in 
(3-|-l) dimension. Investigation on full collapsing model 
of a (2 -f- 1) dimensional star is beyond the scope of this 
analysis. We would like to perform this study in the fu- 
ture. 



FIG. 11: The variation of redshift function Za is shown 
against r at stellar interior. 

VI. CONCLUDING REMARKS 

By taking the Krori and Barua metric as input and 
treating the matter content as anisotropic in nature , 
we have obtained some new types non-singular solutions 
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